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Abstract. Vacuum conditions inside the LTP Gravitational Reference Sensor 
must be under 10 -5 Pa, a rather demanding requirement. The optical window 
(OW) is an interface which seals the vacuum enclosure and, at the same time, lets 
the laser beam go through for interferometric metrology with the test masses. The 
OW is a plane-parallel plate clamped in a Titanium flange, and is considerably 
sensitive to thermal and stress fluctuations. It is critical for the required precision 
measurements, hence its temperature will be carefully monitored in flight. This 
paper reports on the results of a series of OW characterisation laboratory runs, 
intended to study its response to selected thermal signals, as well as their fit 
to numerical models, and the meaning of the latter. We find that a single pole 
ARMA transfer function provides a consistent approximation to the OW response 
to thermal excitations, and derive a relationship with the physical processes taking 
place in the OW. We also show how system noise reduction can be accomplished 
by means of that transfer function. 
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1. Introduction 

LISA Pathfinder (LPF) is an ESA mission, with NASA contributions, whose main 
objective is to put to test critical parts of LISA (Laser Interferometer Space Antenna), 
the first space borne gravitational wave (GW) observatory jl] . The science module on 
board LPF is the LISA Technology Package (LTP) [2] . The unprecedented sensitivity 
of the LTP has prompted the conceptual enhancement of LPF's science objectives as 
regards the purity of geodesic, or free fall motion of test masses in the interplanetary 
gravitational field |S] . 
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Free fall control is achieved by Gravitational Reference Sensors (GRS) [6]. These 
are a set of capacitive sensors which can determine to high precision (nano-metres) 
the 3D position and orientation of cubic test masses relative to their non-contacting 
enclosure, which is rigidly linked to the spacecraft structure. Detected off-centre 
deviations trigger action by a set of micro-thrusters which move the spacecraft such 
that the mass returns to its centred position. The combination of the GRS, the 
thrusters and the control system (the DFACS, Drag Free and Attitude Control System) 
is called drag-free subsystem. The latter is intended to accurately nullify the effects 
of any non- gravitational forces acting on the spacecraft. This makes possible to 
detect differential gravitational accelerations between the two test masses, whether 
by precision interfcrometry [7] or by the drag free system itself. 

This is fundamental for LISA, since Gravitational Waves (GWs) show up as 
tides, i.e., time varying differential gravitational accelerations. The precision of the 
measurement done with the LTP is required to be [5] 



We shall refer to the above frequency band as the LTP Measuring Bandwidth 
(MBW) in the sequel. Equation (JlJ is ten times less demanding than what is 
needed for LISA [5], both in magnitude and in frequency band, yet it is between 
two and three orders of magnitude better than has been achieved or required so far 
for space missions [3]. It has relevant consequences for future missions, which need 
high performance drag free, hence the relevance of LPF beyond its natural objectives 
as a LISA precursor. 

In order to meet the above requirements, the residual pressure inside the GRS 
must be under 10~ 5 Pa, a condition which is classified as very high vacuum by 
the American Vacuum Society. This implies that their interior have to be tightly 




lmHz < — < 30mHz (1) 
2tt 




Figure 1. Layout of the LTP Core Assembly, as of 13-October-2006. The test 
masses are Au-Pt alloy cubes inside either vacuum enclosure. 
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sealed within a vacuum enclosure (VE), and non- mechanical getter pumps installed 
to ensure a suitably rarefied environment around the test masses. Perhaps a more 
obvious option would have been to communicate the VE directly with the external 
interplanetary vacuum, which is much better than 1CP 5 Pa. Recent studies [3] by 
the Project Engineering Team have shown that there are serious difficulties with such 
option. For example, venting out of residual gas has time scales exceeding the very 
LPF mission lifetime, cleanliness control inside the VE is tighter with the OW, etc. 

The general layout is shown in Figure [T] As regards the issues we address in this 
paper, attention is drawn to the Optical Window (OW), which is the interface between 
the test masses and the Optical Bench: laser beams must bounce off the test masses 
to monitor their positions by precision interferometry, hence a transparent window is 
necessary for the light to make it to the interior of the VE. 

The OW is a plane-parallel plate and is therefore a potential source of noise: 
random variations of its optical properties may result in corresponding optical path 
fluctuations, which distort the laser light phase, hence the optical Metrology readout. 
Great care must be taken when manufacturing this critical component of the LTP 
and, once manufactured, characterisation of its behaviour duly performed. 

The most important agent responsible for OW fluctuations is temperature 
fluctuations. These cause various degrees of mechanical stresses across the rim, as 
well as temperature dependent index of refraction changes. The former are very 
difficult to model with quantitative accuracy, mostly due to lack of precise control of 
mounting interface behaviours, but the former can be much better studied in a stress 
free environment. This paper is concerned with the experimental characterisation on 
ground of prototype OWs, and with the phenomenological modeling of their response 
to thermal excitations. This is justified if the noise fluctuations are smaller than the 
applied stimuli and if the system behaves linearly. 

The philosophy of the approach is the one typical of the diagnostics subsystem, 
as described in [10]. This is: apply controlled temperature signals of high signal-to- 
noise ratio to the Titanium flange where the OW is held — see below — , and measure 
the temperature of the former. Measure also the induced phase shifts in a laser beam 
which travels through the OW, then try to establish the transfer function between both 
magnitudes, i.e., temperature and phase-shifts. The transfer function thus obtained 
is assumed to also be valid in the situation when only noise is present in the flange. 
The latter extrapolation hypothesis is the clue to the determination of the phase noise 
contributed by temperature fluctuation noise in the OW, on the basis of temperature 
measurements . 

The success of the proposed empirical approach depends on our ability to find a 
transfer function which depends on a (preferably reduced) number of parameters, 
which do not change significantly across different conditions and runs of the 
experiment. As we shall now show, we have found that a single pole A_RiVL4[|] process 
describes rather satisfactorily the relationship we look for. The precise meaning of this 
concept will be discussed in detail in the following sections, and the results applied to 
evaluate the temperature fluctuation noise in a dedicated experiment. 

The paper is organised as follows: in section [2] we describe the experiment layout, 
including hardware and data acquisition details. Section [3] is devoted to the data 
processing, analysis, two modes of model fitting — a direct linear regression and a 

| Acronym for Auto Regressive Moving Average, roughly the discrete time series equivalent of a 
linear differential equation — see a standard textbook, e.g., [14] 
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Figure 2. The OW (left), with the plane-parallel plate inside the Titanium 
flange, heaters on the lateral surface of the latter (pale brown foils) and two pairs 
of NTC temperature sensors (black beads with wires). On the right, mounting of 
the OW on rails for measurements. Note heaters are covered with aluminium foils: 
this is to prevent undesired heating of other components by heaters' emission of 
thermal radiation. 

single pole ARMA model and numerical results. In section [4] we examine in detail 
the ARMA(2,l) fit, and derive important implications for the understanding of the 
physical processes happening in the OW. Section [5] addresses how the previous analysis 
can be applied to quantify the contribution of temperature fluctuation noise to the 
total phascmeter noise, based on another set of experimental results, and section [6] 
develops an interesting exercise whereby a continuous time model is suggested as the 
origin of the discrete time ARMA fit. Finally, Conclusions and bibliographic references 
close the article. 



2. Experiment description 

The current baseline of LISA PathFindfer and LISA includes vacuum tanks containing 
the test masses which act as end mirrors for the interferometer. Presence of such tanks, 
or vacuum enclosures (VE), force the inclusion of a transmissivc clement interfacing 
between the interior of the VE and the optical bench outside. This optical clement is 
the Optical Window (OW). In this section we describe the laboratory hardware and 
conditions of several runs of measurements conducted in AEI Hannover Laboratory 
facilities to characterise the thermal behaviour of the OW. 

In the experiment two different prototype OWs were tested. Both were 
manufactured following the same baseline as the one to be applied in the final LTP 
flight model. The main element of the window is a very low thermal expansion 
coefficient glass chosen in order to minimise the variation of the optical pathlength 



with respect the temperature. The figure of merit — quantified by equation (211 
below— is of 0.59 x lO^K" 1 for our particular choice, the OHARA S-PHM52 
(n= 1.606, dn/dT— — 5.54 x 10~ 6 K _1 ). This parameter can reach values as high 
as 5.31 x lO^K- 1 for BK7 or 8.32 x lO^K" 1 for fused silica. The glass, of 30 mm 
diameter and 6 mm length was clamped between two Titanium flanges, fastened by 
means of Titanium bolts, and sealed by two helicoflex rings[|] to prevent gas leakage 
in space conditions. 



§ These are softer metallic rings, e.g. Aluminium or Silver. 
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Figure 3. Schematic of the interferometric measurement concept. 

The OW is expected to induce thermal related noise in the Metrology subsystem. 
In order to quantify its contribution to the total noise budget a set of thermal 
diagnostics items were attached to the optical window prototypes. They are shown 
in figure [5] left panel: two Kapton heaters Minco HK5303 attached to the Titanium 
flange lateral surface, and four glass encapsulated thermistors Betatherm G10K4D853 
attached in pairs to the Titanium flange and on the athermal glass surface, for precision 
temperature measurements. These diagnostics items were all glued to their attachment 
points with Pressure Sensitive Adhesive (PSA) tape 3M-966, of similar characteristics 
to the one to be used in flight. The temperature sensors on the glass will actually 
not fly with the LTP. They will however provide relevant information to implement 
real mission data analysis procedures and methods, for which only the Titanium 
temperature data will be available. 

During the experiment, the window was leaning vertically on a PVC two-rail 
structure — see figure [2j right panel — , which impeded any high conductivity thermal 
contact with the rest of the hardware. Although not directly affecting the thermo- 
optical interaction studied here, the OW will be part of the VE in the real LTP, thus 
a higher thermal conductance is to be expected, and therefore a faster suppression of 
thermal gradients is foreseen during mission operations. 

The complete set-up (i.e., the glass plus its mounting structure and the just 
mentioned diagnostics items) was inserted as a transmissive element in a dedicated 
optical bench, as seen in figure [3] The heaters were covered with aluminium foil to 
reduce thermal radiation effects (figure [2] right). For the same reason, the window 
was introduced in a copper box leaving only a narrow opening for the laser beam 
to go through. As seen in the schematic of figure [3] the beam traverses the OW 
only once. This will not be the case in the real LTP, where the laser will go twice 
through each window, instead, but the one passage configuration used here simplifies 
the OW thermal characterisation without information losses. All the experiments were 
performed under low pressure conditions at a P ~ 10 _3 Pa vacuum level. 

The optical window was subjected to various heat pulses comprising a wide 
range of duration length and powers in order to identify suitable parameters for the 
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thermal test to be performed in-flight. The data here reported gather 25 experiment 
runs on two different prototypes, applying heat pulses from 100 mW to 2W ranging 
from 10 s to 100 s of application time. All experiments were performed at room 
temperature, which falls within the expected range of working temperatures of the 
LTP experiment during operations, required to be between 10 °C and 30 °C. Figure [3] 
shows a typical response data plot, with indication of the temperature sensor readings 
and the interferometrically registered phase shifts corresponding to a specific heat 
signal input — see the figure caption for the details. 

Two different data acquisition systems were used in the experiment: the 
intcrferometric data were acquired via the LTP phasemeter prototype [7J, whereas 
the thermal diagnostics data were acquired using the LTP front end electronics (FEE) 
prototype [TTJ [12] . Both acquisition systems have previously successfully passed tests 
of compliance with mission noise budgets. 



3. Fitting the data to models 



The main purpose of this section is to give account of the measured interferometer 
output data in terms of the also measured temperature data. While in this experiment 
both are of course ultimately caused by the heaters' signal, our interest focuses on 
the temperature vs. phase relationship, as this is the one we need to quantify the 
magnitude of temperature fluctuations noise during science operations in flight |10j . 

To serve this purpose, we adopt model fitting techniques. Two approaches will 
be proposed, and discussed in the ensuing section: a direct linear regression fit of the 
intcrferometric data to the temperature read-out coming from sensors on the Titanium 
flange and those on the OW glass itself, and an ARMA model using only temperature 
readout from the Titanium temperature sensors. The latter is of particular interest, 
since it is not foreseen that temperature sensors be attached to the glass surface in 
the real LTP. 




Figure 4. Phase (black curve) and various temperature sensors' responses 
(coloured curves) to a 2 Watt heat pulse applied during 100 seconds. Legend 
indications correspond to the thermometers shown in Figure [2] 
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3.1. Data conditioning 

Before we attempt to fit the data to a useful model, some data pre-processing is 
required. 

The temperature and phase acquisition data systems reside on different hardware 
and software, and deliver the respective time series data for analysis at sample rates 
which arc different as well: temperature data are sampled at / s ,t = 0.65 Hz, whereas 
phase data are sampled at f s s = 32.4 Hz, instead. Downsampling and resampling 
thus needs to be applied to the latter in order to make meaningful sense of data 
fitting algorithms. To avoid aliasing effects at downsampling phase, suitable low pass 
filters are applied. This is however not enough to have matched sampling times in 
both time series, so an additional interpolation algorithm is used for properly matched 
resampling. 

In addition, each data segment is de-trended prior to model fitting. The removed 
trend is evaluated from the first 500 seconds previous to the heat input signal begins. 
This way we get rid of environmental drifts and spurious trending effects. 



3.2. Direct Linear Regression 

A typical phasemeter response when heat pulses are applied to the OW is shown 
in Figure [4j An essentially instantaneous phase response is observed in coincidence 
with thermometers' excitations, which suggests phase behaviour can be described as 
a direct, or single-time relationship between the various temperature readings and 
associated phase shifts. If we additionally make the hypothesis that such relationship 
is Zmear[jj]then the model is given by 

<j>{t) =piT Ti (t)+ P2 T Glass (t) (2) 

where Txi(t) is the temperature read by the thermometer on the Titanium flange 
closest to the activated heater, and TciassW that of a thermometer on the OW glass. 
The parameters pi and p 2 are to be estimated by a least squares algorithm, which 
requires the square error 

N 

e 2 = {<j>[n] - pi T Ti [n] - p 2 Tbi ass [n]} 2 (3) 

to be the smallest possible for the given data streams. Here, 4>[n] is the n-th sample of 
the measured phase, and Txi[n] and Tbiasst^] the corresponding temperature samples. 
Thus, for example, 

<P[n] = <p{nAt) , Tn[n] = T Ti (nAt) , (4) 

etc., where the sampling time At has been set to At = l// s ,Temp, as discussed in 
section [3~T] above. Finally, N is the number of analysed samples of each read-out. 
The conditions of minimum square error are of course given by the two equations 

^— = = ( 5 ) 
dpi dp 2 

which, once solved, give least squares estimates p\ and p 2 of the parameters p\ and 
p 2 , respectively. An example of this procedure is shown in figure ([5]). We report on 



the results of this analysis in section 3.4 



This is in fact quite accurate, on account of the rather small temperature and phase variation 
ranges detected in the experiment. 
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Figure 5. Example of fit results for two different models. 



3.3. ARMA model fit 

Although the linear regression method shows to perform quite acceptably well, there 
is a clear motivation to find a model able to fit the data independently of the glass 
temperature readings, since the latter will not be available in flight. 

In this section we take a different approach to fit phase data to Titanium only 
temperature readings. We shall now assume that the relationship between these 
magnitudes can be expressed by a dynamical equation. More specifically, we make the 
hypothesis that phase is the output of a linear ARMA process^]] whose input is the 
temperature of the Titanium, as recorded by the thermometer next to the activated 
heater. We express this by the formula [15 

<f>[n] = G(q,0)T Ti [n] (6) 

where G(q, 6) is a rational expression of the type 

with q representing the shift operator: 

q x[n] — x[n + 1] , q^ 1 x[n] — x[n — 1] (8) 

for any discrete series x[n\. Finally, 6 is an abbreviation for the vector of r + s+1 
ARMA parameters ao, . . . , a r , /3i, . . . , j3 s . 

System identification in this approach is again based on a least squares criterion, 
for which a suitably defined square error needs to be defined. Following [15j . this is 
the so called prediction error: 

N 

e*(9) = J2Mn]-G(q,e)T Ti ln}} 2 (9) 

n=l 



1 We feel again justified in assuming a linear relationship by the small variation intervals of the 
magnitudes involved. 
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Table 1. Mean values and rms variances of parameter estimates 



DLR 


ARMA(2,1) 


Pi = (-38 ±4) x 10 


~ 3 rad/K 


a = 


(39.6 ±3) x 10~ 3 


rad/K 


p 2 = (65 ±20) x 10" 


3 rad/K 


OL\ = 


(-39.5 ±3) x 10" 


3 rad/K 






01 = 


-0.996 ± 0.001 





The estimates 9 of the parameters 9 are those which cause e 2 (9) to be minimum. 
Algorithms to find them are more robust if the additional hypothesis holds that the 
residuals {(j)[n] — G(q, 9) Txi[n]}, where <j>[n] and Txi[n] are the actually recorded data, 
are a white noise sequence |15j . Reassuringly, this is quite accurately true for our data. 
An example result of the fit is shown in figure [5] blue curve. 

3.4- Numerical results 

As stated in section [2] up to 25 rounds of measurements were carried through 
during the experiment. This is a considerable number which enables us to check the 
consistency of the fitting models just described. The methodology we have adopted is 
the following: for each run, we de-trend the data and then fit them to both the Direct 
Linear Regression (DLR) and the ARMA models. Parameter estimates are then filed 
for further analysis, as we now describe. An observation on the ARMA fit is however 
in order before we proceed. 

Indeed, in the ARMA fit we also need to make a choice of order of the process, 
i.e., we need to set the number of a's and 0's in equation It turns out that 

an ARMA(2,l), i.e., two a's and one (3, is an excellent approach, in the sense that 
differences between model predictions and actual phase data are kept small to a rather 
satisfactory level. Finer tuning can be accomplished adding more a's and/or /3's, but 
only at the expense of excessive parameter estimates' dispersion across different runs. 
This is highly undesirable, and hence we confine our model to the ARMA(2,1). 

Table [l] summarises the results of the analysis, and figure [6] shows the binned 
distribution of the parameter estimates across the 25 runs. An outstanding charac- 
teristic of the fit is that the two MA coefficients very accurately verify the numerical 
relationship a\ = — a®. The model thus relates the output phase data to the time 
derivative of the Titanium temperature — we come back to this in section [4j 

There are a few important aspects of these results which are worth stressing: 

• ai ~ — ao, although the difference between their values is much less than their 
variances. 

• Well within tolerance, —ct\ = a =pi. 

• 0i is strongly peaked at a nominal value, with only 0.1 % relative tolerance. 

• pi is much more disperse, with almost 30% variability. 

4. The ARMA(2,1) transfer function 

In view of the above remarks, it is expedient to rewrite the ARMA(2,1) model in 
terms of the following parameters: 

a = —ai , S = ao + a.\ , = 0\ (10) 
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is the z-transform of the process transfer function — we have replaced the shift operator 
q by the complex variable z [TJ. It is also expedient to emphasise the structure of 
this formula as follows: 

G(z, a, 13,5) = a G HP (z, (3) + 5 G L p(z, (3) 

with 

1-z- 1 „ , 1 



G UP (z,[3)= 1 , _ 



, , . | G LP (z,P) = 
1 + pz 1 1 + pz 1 



We now find discrete Fourier transforms (DFT) by the substitution 
z = exp(iuAi) 

where At is the sampling time of the time series data. The following obtains: 

4sin 2 (wAi/2) 



|G HP (w,/3)| 2 
\G LP (^,(3)\ 2 
\G{u,a,P,5)\ 2 



1 + 2/3 cos(wAi) + /? 2 
1 

1 + 2/3 cos(wAi) + (3 2 

+ 5) sin 2 (wAt/2) 
1 + 2/3 cos(wAi) + /3 2 



(12) 
(13) 
(14) 

(15) 
(16) 
(17) 
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a parameter [rad/K] 



20 



-1.002 -1 -0.998 -0.996 -0.994 -0.992 -0.99 
(3 parameter 
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0.05 0.1 J „^ 0.15 

p parameter [rad/K] 



Figure 6. Histograms of the estimated parameters for the two fitting models 
described in the text. Top: ARMA(2,1), with ol = olq and j3 = j3\. Bottom: DLR,. 
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Figure 7. Bode diagram for the optical window ARMA(2,1) transfer function 
estimate using the values from Table]!] Left panel: its modulus; dashed lines show 
the 1 a confidence region. Note that this confidence region has been extrapolated 
below ~0.4 mHz, as actual experimental data were actually not available in that 
band. Right panel: phase of the transfer function. 



We thus see that the transfer function G is the sum of a high pass term, Ghp , and 
a low pass term, Glp. The effect of the latter is naturally dominant at low frequencies, 
while the high pass term dominates at high frequencies. The concepts of low and high 
frequencies can be made precise by means of some intermediate frequency /b where 
the gains of Ghp and Glp are equal. This is easily calculated: 

' (18) 



fb 



2nAt 



and has a numerical value of /b — 0.2 mHz, which means the high pass effect dominates 
throughout the LTP bandwidth. We may not however neglect the relevance of the low 
pass at lower frequencies, as it contributes extremely valuable information for LISA. 

A Bode diagram representation for the transfer^functions is shown in figure [7] 
The filter modulus is characterised by a plateau of |G| ~ 40 x 10 -3 rad/K across the 
entire LTP measuring bandwidth. Temperature fluctuations at frequencies below this 
bandwidth arc also suppressed but following a different behaviour, related to the low 
frequency response of the optical window. The Figure also shows the phase behaviour 
of the filter, and indicates increasing delays for high frequency perturbations. 

The DC gain of the filter is 

\G(u = a,a,[3,8)\ = ^ (19) 

If the estimated parameters are substituted in this expression then the following 
is obtained: 

\G{uj = 0, a, p,S)\ = (25 ± 4) x 10" 3 rad/K (20) 

We now try to produce some insight into the physical meaning of the just discussed 
analysis facts. 



4-1. Physics of the ARM A process 

Two different kinds of thermal effects have been identified as sources of changes in the 
optical path-length of a light beam traversing a plane-parallel piece of glass: 
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i. Temperature dependent changes of the refraction index 

ii. Mechanical stress induced changes of the refraction index 

We briefly describe below how these effects can be approximately evaluated. 



4-1.1. Temperature dependent changes of the refraction index The first effect, which 
is found under stress free conditions, is quantified by the formula |16j 



d<t> 
dT 



free 



L 

= 2 *A 



dn 
dT 



+ (n — 1) ole 



(21) 



where <\> is the phase shift suffered by a beam of light traversing a glass slab of thickness 
L and (nominal) index of refraction n; A is the wavelength of the used light, and cue is 
the linear thermal expansion factor of the glass, ole — L dL/dT. 

The d(j)/dT\f rcc effect is most prominent at very low frequencies and DC. The 
reason is that it happens even if the temperature of the glass is homogeneous, and 
without mechanical stresses. It has been measured on naked glass samples in the 
laboratory, free of any pressure or tension, with the result that it is 25mrad/K [17], a 
figure very well matching the one given by equation ( 20 1 . One however needs to 
consider that the latter was obtained from data of a real window, i.e., including 
metal flange. This consequently means that the stress contribution dip /da must be 
comparatively small at very low frequencies. 

The same result is endorsed by another independent evidence. If data-sheet 
properties of the OHARA S-PHM52 glass used in the experiment are used to calculate 
the thermal related path-length variations in the optical window glass due to changes 
in the refractive index, the result is that d(j)/dT\{ iee is ~21mrad/K, again in good 
agreement with equation ( 20 1 . 



4-1-2. Mechanical stress induced changes of the refraction index This second effect 
is relevant to our experiment because the glass, clamped by Titanium flanges to the 
ISH structure, is under stress due to differing thermal expansion coefficients in glass 
and metal. Mechanical stress also induces pathlcngth changes which are difficult to 
model. From the datasheet, the only parameter provided by the manufacturer which 
can be used to quantify these interactions is the photoelastic coefficient, (3. However, 
it must be noticed that (3 does not describe the change in the refraction index due to 
stress, dn/da, but the appearence of birefringence due to stress, i.e., the change of the 
velocity of light along different axes of the material. Although not directly related, 
both parameters range in the same order of magnitude [TB], and we shall thus use 
the photoelastic coefficient here for our order of magnitude estimate, described in the 
following. 

Under this simplyfing assumption, the photoelastic coefficient can be related to 
a pathlength variation by 

Assess = /3ad (22) 

where (3 = 10~ 5 nmcm -1 Pa -1 , d is the glass thickness (d = 0.6cm for the Optical 
Window), and a is the applied stress, having dimensions of pressure. 

In this case, the stress on the glass is due to differential thermal dilatation of the 
Titanium flange and the OW glass itself. The situation is illustrated graphically in 
Figure [8] Because the coefficient of thermal expansion of the glass is larger than that 
of the Titanium flange embracing it, the latter expands less when submitted to the 
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Titanium alone 




Glass 
alone 



equilibrium 

Figure 8. Schematics of the dilatation of the OW glass and the clamping titanium 
flange: the right part zooms in the profiles acquired by the interface (red) when 
the temperature of the block increases by AT. Should either the Glass or the 
Titanium be let expand freely, the boundaries would be placed as also represented. 
Dashed line is the interface position before heating. 



same temperature rise, and hence the glass is compressed radially along the rim. The 
opposite happens if the temperature decreases, i.e., the glass is in this case stretched 
outwards by the radial pull of the Titanium. The contraction/expansion forces acting 
on glass and Titanium reach an equilibrium state which determines the radii of the 
contracted/expanded pieces of Titanium and Glass. The equilibrium position thus 
happens when 

[Sp T + <W] xi = i S PT - fy<r] Glass ( 23 ) 

where 6pr and 5p a refers to changes in radius caused by temperature changes and 
by stresses, respectively. The above formula holds even if temperature changes in 
Titanium and glass are unequal. On the other hand, we are not considering in our 
description possible effects coming from the helicoflex ring between the Titanium and 
the glass. As stated above, we are here trying to get an order of magnitude of the effect 
based on a simplified mechanical model, and interface effects are thus not included. 



The contributions appearing in equation ( 23 ) are given by |18j 
SpT = paAT and 5p a 



VP 
IE 



(24) 



where p is the radius of the interface between Titanium and the glass, £ stands for the 
width of the body, E is the Young modulus, a the thermal expansion coefficient and 
p the lateral pressure. Combining equations (23 1 and (24 1, and following the notation 



of figure [8j we find the lateral pressure on the glass: 



P : 



axiATxi — aciassATciass 



Er 



(r/h)E Ti -r^ Glass 

The strain on the glass lateral surface is given by o-Qia 
is the thickness of the window glass — see [TH]. Hence, 



(25) 

= prdj [t d) = p, where d 

(h/r)E~l ss (26) 

We can consider two different regimes here: the low frequency (LF) regime and 
the high frequency (HF) regime. The first corresponds to long duration heat pulses 



_ QTiA.r Ti — a G iassAr Glass 

O'Glass — 1 , ,, , N „_i 

Ti 
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applied on the Titanium flange, actually long enough that the temperatures of the 
glass and the Titanium equal each other, or ATx; = ATciass = AT. In this case 

axi — aciass 



AT , low frequency 



(27) 



E^ + (h/r)E G l 

On the other hand, if short heat pulses are applied on the Titanium then the glass 
does not have time to respond, and in this we can assume ATciass = and ATti = AT. 
Thus, 

C*Ti 



OGla 



E^ + (h/r)E G l 



AT , high frequency 



(28) 



We can use the above formulas in combination with (22 1 to obtain 

2nd 



dT 



Stress 



B 



B 



laser 

2nd 



axi — QfQlass 

(h/r)E G l 
■ i i 



E, 



Ti 



Glass 



laser 



E- 



Wr)E G l 



low frequency 



high frequency 



(29) 



J Ti "T" VI ' I -^Glass 

It is recalled that A(j> = 27rAs/Ai asor , where Ai aser is the laser wavelength. We put 
numbers here: 



1(T 3 x 1CT 9 Pa" 1 

-3 



d = 6 x 1CT " m 
Ai ascr = 1.064 x 10- 6 m 

a Ti = 8.6 x 10- 6 K- 1 
E Ti = 11.6 x 10 10 N m- 2 
h = 0.02 m 

a Glass = 10 x 10- 6 K- 1 
^Giass = 7.15 x 10 10 N m- 2 
r = 0.015 m 



to obtain 



d4_ 
dT 



Stress 



2.5 x 10~ 3 rad K" 
15 x 10^3 rad K" 



low frequency 
high frequency 



(30) 



4-.1.3. Discussion of the results The total thermal effect is the sum of the above two 
effects, i.e., optical pathlength changes induced by pure thermal expansion and by 
mechanical stress. The former gives a value of 21 x 10~ 3 rad K _1 throughout the 
frequency band, as extracted from datasheet values — see section |4.1.1| above. We 
can thus summarise the results as shown in Table [2] 



Table 2. drp/dT, units in mrad K -1 . 

ARMA Analytic 

LF range 25 ± 4 23J5 

HF range 40 36 

Nude glass 21 



The agreement between the results produced by our simplified model and the 
ARMA fit is quite good. Even though the model is not fully comprehensive of all the 
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physical effects happening in the OW, it can be considered rather satisfactory from 
a purely empirical point of view, hence very useful for practical purposes. Work is 
currently in progress for a more thorough approach, and we shall report on new results 
in due course. 

We conclude from this discussion that the low pass component of the transfer 
function is almost exlusively related to the d(f>/dT\i ICC effect, while the stress effects 
only show up significantly in the higher frequency band. This makes sense, as stresses 
applied along the glass rim quickly propagate inwards throughout the glass piece. 

Although the LTP spectrum is only above 1 mHz, an analysis at frequencies below 
this one, down to 1CP 4 Hz and even further, must be considered of high interest, as the 
latter frequency band will be important for LISA. The experimental data reported in 
this paper can be improved to access the lower LISA band, since they typically consist 
in one hour long runs. This is a strong suggestion for the LTP experiment plan. 



5. Noise projection 

One of the main scientific objectives of the diagnostics system in the LTP is to measure 
identified environmental disturbances [5], and to provide the data and analysis tools 
to estimate the contribution of those disturbances to the overall mission noise budget, 
equation ([lj. In practice this means the LTP Data and Diagnostics Subsystem (DDS) 
must be able to provide suitable transfer functions to convert measured disturbance 
noise into test mass acceleration noise. This section is devoted to describe this 
procedure in the case of temperature fluctuation noise in the OW, and to show how 
it works in an on-ground laboratory experiment — to be extrapolated to a space- 
borne one. 

We will use the results derived in the previous analysis to obtain an estimation for 
the thermal contribution to the interferometer performance. We shall naturally limit 
ourselves to ARMA model, since it is the one making sense for real mission purposes, 
as already discussed. 

The basic idea is that the OW transfer function, as determined from high SNR 
system response, also applies when there is only (weaker) noise in the window |10j . 



For this we shall use the one in equation ( 11 ), i.e., 



G ^^ = a YT0 + YTl^ (31) 

We now show which procedures must be applied to address the problem of finding 
the contribution of temperature fluctuations noise in the OW to the total OW noise. 
To this end we consider data of temperature and phase noise generated in a different 
experiment, and apply to it the methodology just sketched. 

The laboratory setup and the experimental details can be found in reference [19] . 
In this experiment, the optical window is not part of a testing optical bench (OB), 
but is glued in a lateral side of the LTP OB Engineering Model, instead. This way, a 
double beam pass across the window is forced: the laser light travels from the optical 
bench through the OW to a dummy mirror, faking a test mass; there, it is reflected, 
sent back again across the OW and out to the OB. Such setup proved to be compliant 
with the interferometer noise budget, showing that the inclusion of the OW does not 
degrade the interferometer performance. Two temperature sensors in the Titanium 
flange and one in the glass were left in place to measure temperature values during 
long term runs. 
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Figure 9. Temperature fluctuations spectrum (left) measured during the 
performance experiment compared to the LTP temperature requirement and the 
lower limit set by the front end electronics sensitivity limit. Phase fluctuations 
(blue curve, right panel) during the same run are compared to the optical window 
thermal contribution (magenta curve, right panel), as derived with the ARMA 
transfer functions given by equations | |32| l and \17\ . The la confidence region is 
also included for the latter. 



No thermal disturbances were deliberately introduced, so the thermometers only 
read environmental temperature fluctuations. We use equation (17 1 to convert the 
temperature flu 
Sl'fa). Thus, 



1 /2 

temperature fluctuations spectral density, S T (a;), into a phasemeter spectral density, 



Sjg(w) = 2 \G(u,a,0,S)\ 4 /a (w) (32) 

where the numerical factor 2 is required to account for the double passage of the laser 
beam through the window in this case. We assume both passages are completely 
correlated, given the extremely small time scale of their occurrence compared to 
thermal reaction times. Spectral densities are therefore linearly added. 

Results obtained in a typical run with that setup are plotted in figure [9] using 
a MATLAB package being developed ex professo for the LTP Data Analysis [20 . 
The left panel shows temperature fluctuations measured in the Titanium flange. 
As can be seen, these reach the front end electronics (FEE) floor noise in the 
higher frequency region of the measuring bandwidth, while keeping slightly above the 
LTP maximum temperature fluctuations requirements limit in the lower frequencies, 
around 1 mHz |21j . This is in fact a worst case condition, since the temperature 
power spectrum decreases as frequency increases, and thus if the LTP temperature 
requirement is reached at the lower frequency range then the rest of the spectrum will 
naturally follow a descending curve like the one shown in figure [9j 

The phasemeter fluctuations spectrum is however below the required noise level, 
as we can see in the blue curve of the right panel. The temperature fluctuations 
spectral data in the left panel arc then submitted to the algorithm, equation (32), and 
the result is the magenta curve displayed in the right panel. 

The low coupling to thermal disturbances implied by G(uj, a, (3, 5) causes the 
thermal contribution to only represent 5 % of the phasemeter noise at 1 mHz, and 
about 0.5% of the LTP goal. We thus feel reassured that there is still considerable 
margin here. 
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6. Continuous time models 

On account of the empirical results reported in section [3j and of the remarkable 
accuracy with which both DLR and ARMA models fit the experimental data 
— notwithstanding their completely different nature — , we now try to shed some light 
into the kind of processes which take place in the system. 

For this, we attempt to picture the ARMA(2,1) model relating the phase readout 
of the interferometer <p[n] and the temperature at the Titanium flange Tr^n] as the 
digital implementation of some analog physical process. The starting point is of course 
the digital algorithm, equation Q, which in this case is given by 

G(q, a, 0,5) = a + TT^pi ( 33 ) 

where q is the shift operator of equation Q . The recursive form of the process thus 
defined is therefore 

<t>[n] + P4>[n - 1] = a{T Ti [n] - T Ti [n - 1]} + ST Ti [n] (34) 

and can be regarded as the digital implementation of a first order continuous time 
filter, governed by a first order differential equation: 

j>(t) + t- 1 0(f) = Af Ti (t) + BT Ti (t) , f= j^j (35) 

where r is the characteristic time constant of the analog filter, and A and B are scale 
factors, respectively weighing the contributions of the temperature's time variation 
rate and the temperature itself to the phase shift effect. We have assumed the Txi(f) 



dependence in the rhs of ( 35 1 in line with the fit result expressed by the rhs of (34 1 



If the time constant r is much larger than the sampling time At implicit in 



equation (34 1 then we can approximate time derivatives by 



m * «"-«'-*> (36) 

and, mutatis mutandi, the same for Txi(f). Taking t = nAt for the timing of the n-th 



sample, and using the natural notation <fi[n] = 4>(nAt), equation (351 is approxi- 
mated by 

4>[n]-(l + ^\ 4>[n-l]=A(l + ^j {T Ti [n]-T Ti [n-l]} 

+ BAt(l + —) T Ti [n] (37) 



This can be readily compared to equation ( 34 1 to obtain 

/9 = -( 1 + ^)-\ a = A(l + ^y\ 5 = BAt(l + ^-y 1 (38) 

P is seen to have a value very close to —1 (Tablc[T|), or j3= —(1 — rj) with rj < 1CP 2 
comfortably in all cases. Hence r~ At/r], i.e., r> At, which a posteriori justifies the 
approximation leading to equation (|37[) 




The formal solution to equation ([35]) can be easily written down. After initial 
transients die out, the phase is given by 

ct>{t)=AT Ti {t) + (B-A)T- 1 [ e-< t -*'>/ T T T i(0# / (39) 
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The meaning of this filter equation is better understood if we recast it in frequency 
domain: 



H = 



A 



1 



B 



1 



IUJT 



T T1 (w) 



(40) 



This equation shows again that the analog process is also the superposition of 
two contributions: a high-pass filter proportional to A, and a low-pass contribution 



proportional to B: the first arises in equation (34 1 due to the Titanium temperature 



derivative, while the second appears related to the term proportional to the Titanium 
absolute temperature. This split dependence of the OW response to temperature 
pulses points to two different physical thermal processes affecting the glass, as already 
discussed in section SHI 



We can now make use of equations (38) to identify the coefficients A and B in 
terms of the fit parameter values of Tablejl Taking Ai/r<C 1, we find that A ~ a, 
and B ~ 5/ At. In addition, we can take advantage of the relationship a ~ p\ between 
the auto-regressive and the DLR model parameters to obtain an expression relating 
both models. Accordingly, equation (39 1 can be rewritten as 

<j>(t) ~ p x T Ti {t) + (6 /At - pi) r- 1 f e -(*-*')/r Tt .^) dt > (41) 

Jo 

If we go back to the DLR fit formula, equation the following expression 
ensues: 

rt 

e 

la 

after the term 6/ At has been been safely neglected in front of p\. We thus see that 
temperatures in the Titanium flange and in the OW glass are related by a low-pass 
with a time constant, r, of a few hundred seconds — note that p\ and pi have different 
signs, Table [I] 

It must be recalled that this relationship emerges out of the good quality of the 
fits by both DLR and ARMA(2,1) models, and is key to understanding why only 
the Titanium gauge is required to make a good prediction of the OW response to 
temperature variations, as will be required in flight. The physical reason for the 
observed relationship between temperatures is to be sought in the properties of the 
interface between the Titanium and the glass in the OW. 



Iblassf*) — — — T 1 
P2 



r T .( t ') dt' 



(42) 



7. Conclusions 



While the optical window is a crucial element in the LTP optical Metrology system, it 
thankfully appears that it is quite stable to temperature fluctuation noise — so far as 
the latter is compliant with mission environmental requirements. The present paper 
contains a rather thorough analysis of such behaviour, based on experimental data 
gathered through different runs of on-ground laboratory measurements. 

Our main purpose was to prepare for thermal diagnostics analysis tools in 
flight, and to gain as much understanding of the underlying physical processes as 
possible. This means we need to know how noisy data retrieved by thermometers 
can be converted into phasemeter fluctuations, thereby quantifying the contribution 
of temperature random variations to the total mission noise budget — which is the 
ultimate objective of LPF in preparation for LISA. Our most relevant finding is 
the discovery that temperature readings in the Titanium flange embracing the OW 



Thermal diagnostic of the Optical Window. 



19 



plane-parallel plate relate to phase values through an ARMA{2,\) transfer function. 
Although this is the result of numerical analysis, hence lends itself to parameter 
estimation variances, it appears to be considerably robust. 

The analysis has shown that the ARMA(2,1) process naturally splits up as the 
sum of a high-pass and a low-pass process, each of them with significantly different 
relative weights which result in the high-pass dominating above ~1 mHz, while the 
low-pass takes over in the lower LISA frequency band, i.e., at 0.1 mHz and below. A 
major achievement of the analysis has been the identification of the physical processes 
responsible for this behaviour: mechanical stresses — induced by differential thermal 
expansion of metal and glass — are associated to the high pass term, while d(j)/dT\f lee 
effects account for the low-pass. 

We consider the analysis presented here as rather complete in some of its essential 
traits. But there are still open issues which call for further study. For example, heater 
generation of test signals must be monitored by temperature sensors close to the 
activated heaters — due to lag effects in remoter spots — for the procedures described 
herein to be fully operative. This raises some caveats regarding full applicability of the 
noise projection algorithms, as the sources of heat dominating a given temperature 
reading may not be clear in LTP science operation mode. 

A more global tool for full LTP thermal diagnostics, which takes into account the 
specific features of each individual part of the system must be assembled. Research 
on this is currently underway which will reported in due course. 
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